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Homotopy perturbation and Elzaki 
transform for solving Sine-Gorden and 
Klein-Gorden equations 


E. Hesameddini and N. Abdollahy 


Abstract 


In this paper, the homotopy perturbation method (HPM) and Elzaki 
transform is employed to obtain the approximate analytical solution of the 
Sine Gorden and the Klein Gorden equations. The nonlinear terms can be 
handled by the use of homotopy perturbation method. The proposed homo- 
topy perturbation method is applied to reformulate the first and the second 
order initial value problems which leads to the solution in terms of trans- 
formed variable, and the series solution that can be obtained by making use 
of the inverse transformation. 


Keywords: Homotopy-perturbation method; Elzaki transform; Sine-Gorden 
equation; Klein-Gorden equation. 


1 Introduction 


In this paper, we have considered the Sine-Gorden (SG) and Klein-Gorden 
(KG) equations, as the following and also in [2] respectively: 


Utt — Uxa + ag(u) =? f (B20); (1) 


and 
Utt — Ura + Bu + Bog(u) = (ed); (2) 


where u is a function of x, t and g is a nonlinear function . The a pa- 
rameter is so-called dissipative term, which is assumed to be a real number 
with a > 0. When a = 0, Eq. (1), reduces to the undamped SG equation, 
and when a > 0, to the damped one. f is also a known analytic func- 
tion. The Sine-Gorden and Klein-Gorden equations model many problems 
in classical and Quantum mechanics, solitons, and condensed matter physics 
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[1, 5]. The SG equation arose in a strict mathematical context in differential 
geometry in the theory of surfaces of constant curvature [17]. For the SG 
equation, the exact soliton solution have been obtained in [16], using Hirotas 
method in [21], using Lambs method in [18], by Backlund transformation 
and painlevé transcendents in [18]. Numerical solutions for the undamped 
SG equation have been given among others by Guo et al. [8] by use of two 
diffirent schemes, Xin [24] who studied SG equation as an asymptotic reduc- 
tion of the two-level dissipationless Max well-Bloch system, Christiansen and 
Lomdahl [3] used a generalized leapforg method and Argyris et al. found 
the accurate and efficient methods for solving such equations which is an 
active research undertaken by Herbst et al. The method [14], presented a 
numerical solution for the SG equation obtained by means of an explicit 
symplectic behavior of a double-discrete, completely integrable discretization 
of the Sine-Gorden equation, and they have illustrated their technique by 
numerical experiments. Wazwaz [23] has used the tanh method to obtain 
the exact solution of SG equation. Approximate analytical solution of Kline 
Gorden equation thrangh the Adomian Decomposition Method (ADM) was 
presented in [4, 6, 20]. Kaya has applied the modified ADM (MADM) for 
obtaining the approximate analytical solution of the Sine-Gorden equation 
n [19]. Another more powerful and convenient analytical technique, called 
the homotopy-perturbation method (HPM), was first developed by He [13]. 
Some part of this work on HPM can be found in [9, 10, 11]. HPM trans- 
formes a difficult problem into a set of problems which are easier to solve. 
Chowdhury and Hashim [2], have used the HPM to obtain the approximate 
analytical solution of Sine-Gorden and Klein-Gorden equations. Recently, 
Tarig Elzaki [7], has introduced a new integral transform, named the Elzaki 
transform, and it has further applied to the solution of ordinary and partial 
differential equations. 


Now, we consider in this work the effectiveness of the homotopy-perturbation 
Elzaki transform method to obtain the exact and approximate analytical so- 
lution of the Sine-Gorden and the Klein-Gorden equations. 


2 Elzaki Transform 


The basic definition of modified form of Sumudo transform or Elzaki trans- 
form is defiend as follow, Elzaki transform of the function f(t) is: 


EO vf f(t)e~/dt, t>0. (3) 


Tarig M. Elzaki showed the modified form of Sumudu transform or Elzaki 
transform in which it is applied to the differential equations, ordinary dif- 


Homotopy perturbation and Elzaki transform for ... 35 


ferential equations, system of ordinary, partial differential equations, and 
integral equations. 


To obtain the Elzaki transform of partial derivative, we use the integration 
by part, and then we have: 


Of (x,t), 1 
E| ot me ple) arf (e0); (4) 
OP f(z,t),_ 1 Of (x,t) 
E| Ot _ art (zu) f(é,0) =v Ot ’ (5) 
Of(a,t),  d 
ro = ag (&¥); (6) 
Pf(z,t)) _ 
[ Ox2 dx2 T(x, v). (7) 
Proof: First, we assume 
| rae dt = T(x, v). (8) 
By using the integration by parts one obtain: 
Of (a, t) =. Of = nae = Of 
E| yi ]= . Uae eae te . ve att 
—t Dey 
= lim {lve* fle, 18 — f° e+ fle, 1)ab} 
poo 0 
T(2,v 
= TE) _ 72,0), (9) 


Assuming that f is a piecewise continuous function and er exist, also it is 
of exponential order which is means that there exist nonegative constants M 
and T such that for all t > T, we have 


|f(t)| < Me*. 
Now, 
pic i yOFlet) = a6 
Ox 0 Ox 
using the Leibnitz rule to find: 
Of(a,t),  d 
E| au |= ag l(t): (11) 


By using this method, we have: 


pi Plo) a 
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To find: 82 f(x, 1) 
E|—3 2 — |. 13) 
Let: af («,t) 
DE =4%9; 14) 
then we have: 
Pf(z,t),_ ,,O9(@,t), _ _[9(a,t)] 
= B= ee ale,0), (05) 
O* f (x,t iL Of (x,t 
EAD) = Save,v) — f(0,0) — vAS2 (16) 


Therefore, one can easily extend this result to the nth partial derivative by 
using mathematical induction. We will see that, the Elzaki transform rivals 
the Laplace transform in solving the problem, its main advantage is the rivals 
that it may be used to solve problems without resorting to a new frequency 
domain because it preserve scales and units properties, the Elzaki transform 
may be used to solve intricate problems in engineering, mathematics and 
applied science without resorting to a new frequency domain. 


3 Homotopy Perturbation Method: 


The basic idea of the standard HPM was given by He [9, 12] and a new 
interpretation of this technique was presented by our research gorup [15]. 
To introduce HPM, considered the following general nonlinear differential 
equation: 

Lu+Nu= f(z,t), (17) 


with initial conditions: 
u(@,0)= ky, ue (a, 0) = ko. (18) 


Where u is a function of x, t and k1, kp are constants or functions of x. Also 
Land N are the linear and nonlinear operators respectively. According to 
HPM [2] we construct a homotopy which satisfies the following relation: 


H(u,p) = Lu— Ly + p[Luy + Nu — f(a, t)] = 0, (19) 


where p € [0,1] is an embedded parameter and 1 is an arbitrary initial 
approximation satisfing the given initial condition. 
By setting p=0 and p=1 in Eq. (19), one obtain: 


H(u,0) = Lu-—Lv = 0, and H(u,1) = Lu+ Nu-— f(a, t) =0, (20) 
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which are the linear and nonlinear original equations, respectively. In topol- 
ogy, this is called deformation and Lu— Lvg and Lu+ Nu-— f(a,t) are called 
homotopic. Here, the embedded parameter is introduced much more natu- 
rally, unaffected by artificial factor, further it can be considered as a small 
parameter for0<p<1. 

Chowdhury and Hashim in [2], have presented an alternative way for 
choosing the initial approximation, that is: 


Vo = u(x, 0) + tuy(x,0) + Lo" f (x,t) = ky + tke + L' f(z, t), (21) 


where L~!(.) = ‘fe te i Jo ()at...dtdt depends on the order of the linear op- 
erator. With this assumption that the initial approximation vp given in Eq. 
(21), in HPM, the solution of Eq. (19), is expressed as: 


u(x, t) = uo(a, t) + puz (a, t) + p?ug(a,t) +... (22) 


Hence the approximate solution of Eq. (17), can be expressed as a power 
series of p, i. e. 


u= lim w= Sou. (23) 


4 Homotopy Perturbation and Elzaki Transform 
Method: 


Consider a general nonlinear partial differential equation with initial condi- 
tions of the form: 

Du(z,t) + Ru(z,t) + Nu(2,t) = flv, t), (24) 

u(#,0) =c1, ur(x,0) = co. (25) 


Where D is a linear differential operator of order two, R is a linear differential 
operator of less order than D , N is the general nonlinear differential operator, 
f(x, t) is the source term and c1,c2 are constants or functions of x. 

By taking Elzaki transform to both sides of Eq.(24), result in: 


E|[Du(2,t)] + E[Ru(a, t)] + E[Nu(a,t)] = E[f(a,t)). (26) 


Using the differentiation property of Elzaki transform and the initial condi- 
tion in Eq.(24), one obtain: 


Elu(a, t)] = v? Elf (a, t)] + ve, + v'cq — v° E[R(a,t) + Nu(a,t)]. (27) 


Applying the inverse Elzaki transform on both sides of Eq.(27), we get: 
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u(az,t) = F(2,t) — E~\{v? E[Ru(z, t) + Nu(2, t)]}, (28) 


where F(x,t) represent the term arising from the source term and the pre- 
scribed initial conditions. 


According to HPM we have: 


u(x,t) = F(w,t) — pE“{v®B[Ru(a, t) + Nu(z,2)]}, (29) 
now by substituting 


i=0 i=0 
in Eq. ( 29 ) where H;(u) is He’s polynomials that are given by: 
H, hy a 01,2 31 
(Uo, U1, --) Ui) = Op OUP Us) p04 tO Gly 


one obtain: 
i=0 


— PEP ERD p'ui(e,t)) +N) p'Hi(u))]]}- (82) 


1=0 i=0 


This our method is infact a coupling technique of Elzaki transform and the 
homotopy perturbation method. Comparing the coefficients of the like powers 
of p the following approximations are resulted: 


D : uo(x, t) F(a; 2); 


p' : uy(2,t) = —E~'{v? E[Ruo(a, t) + Ho(u)]}, 
p* : Ue(2,t) = —E~'{v° E[Ru, (2, t) + Hi (u)]}, 


I 


I 


etc. 
Therefore the solution will be obtained as: 


u(x,t) = uo(a, t) + ui(a, t) + ue(a,t) +... (33) 
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5 Numerical Applications 


In this section, we apply the homotopy-perturbation and Elzaki transform 
method for solving Sine-Gorden and Klein-Gorden equations. Numerical re- 
sults are very encouraging. 


Example 5.1 First, consider the following Sine-Gorden equation with the 
given initial conditions: 


Ut — Ure +sinu=0, u(av,0)=0, uz(x,0) = 4sech(x). (34) 
The exact solution is given in [2] as: 
u(a,t) = 4arctan|[tsech(x)]. (35) 


= 
ue 


To solve the example by this method, we take sinu ~ u — - Ata: 


After taking Elzaki transform of (34), subjected to the initial conditions, 
one obtain: 


E[(u(z, t)] = 4u°sech(x) + Blut, t)] — v?Efu | 


The inverse Elzaki transform implies that: 


u(x, t) = dtsech(x) + Be? “SBI (u(x,t)] v°Efu EG ale (37) 


Now applying the homotopy perturbation method, we get: 


co 2 co 
Do viulost) = ascents) +(e? PID pints) (38) 


(CEopiule,t)® | OEoriu(e.t)® 


6 120 I}. 


3 VED, p'u,(a, t) 
i=0 


By comparing the coefficients of the same powers of p, result in: 


p® : ug(a, t) = 4tsech(x), 
p':us(a,t) = © 1, 44sech(x)u?E(t)} — E~'{4tsech(x)v?E(t)} 


dx? 
1024 


sf BM sech® (x)? E(t%)} — B-1{ = sech®(x)vE(t*)}. 


Then, we get: 
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p' : u,(a,t) = —~sech*(x) — (—105t? cosh?(x) + 42t® cosh?(x) — 16t”), 
4 Boach? 8 9 2 4 4 
: = ———_— 40t° — t h — 42 h 
pr : U2(a, t) 5007005" sech” (x)(7040t° — 33696t” cosh“ (x) 90t* cosh” (x) 
+ 143000¢* cosh? (a) — 205920t? cosh* (a) + 514802? cosh® (a) 
— 270270 cosh®(«) + 405405 cosh* (ar), 


therefore, the 3-terms Elzaki-HPM solution is: 


4 
u(x,t) = s0n7005 sec (x)(7040t!? — 33696t?° cosh?(x) — 4290t® cosh*(x) 


+ 143000¢° cosh?(«) — 308880¢° cosh* (a) + 51480t° cosh® (2) 
+ 405405t* cosh*(a) — 675675t? cosh® (x) — 270270 cosh® (x). 


The behavior of the solution (34), by 4-terms of EHPM and its exact solution 
are shown in Figures 1 and 2. 


Exact solution 


Figure 1: These surfaces show the approximate solutions obtained by 4-terms 
of EHPM and the exact solution of u(x,t), respectively .(a) Exact plot; (b) 
EHPM plot(Eq. (35)) 


Example 5.2 Considering the following Sine-Gorden equation with the given 
initial conditions: 


Utt — Ure + sinu = 0, u(z,0)=a+ecos(uxr), u(av,0)=0. (39) 


Where pp = a and ¢ is a constant. 
We take sinu ~ u — = + — to solve this example. 
By taking Elzaki transform of (39), subjected to the initial conditions,we 


have: 


2 3 5 


E|u(z, t)] = v?(m + € cos(ux)) + v? _ E[u(z, t)] — v? Elu 
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x x 
E-HPM E-HPM 


(a) (b) 


s+ * Exact 


crt Exact 


Figure 2: Comparison of the results by EHPM with the exact solution for 
example 5.1. 


Table 1: Comparison of the exact and numerical solutions of example 5.1 


x t Exact solution EHPM-4term Absolute error 
0 0 0 0 0 

0.062 0.162 0.6412071928 0.6412071885 4.3 x 1079 
0.156 0.156 0.6116706532 0.6116706503 2.9 x 107° 
0.469 0.469 1.5964486815 1.596478333 .8482 x 1077 
0.539 0.403 1.349539812 1.349538454  .1358 x 107§ 
0.620 0.685 2.077035544 2.076989415  .0590 x 1079 
0.781 0.781 2.136014353 2.135871203 .1431 x 1077 
0.896 0.893 2.234155992 2.233473035  .6829 x 1077 
0.975 0.992 2.319920856 2.317982103  .1938 x 1077 
1.000 1.000 2.300024730 2.297874130  .2150 x 1077 


The inverse Elzaki transform implies that: 


u(a, t) = (w +e cos(ux))E~*[v?] + Bu? SS Bl(u(a, t) —v Elul 


3 
fy 


2 [u*] 


To! 


Now applying the homotopy perturbation method, result in: 


u(x,t) = (7 + €cos(px)) EB *[v?] + p{EBMv? SS Bl(u(a, 1) —v E[ul 


3 
2 Uu 
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By taking u(x,t) = S072 p'ui(x,t) and comparing the coefficients of the 
same powers of p, one obtain: 


p° : uo(z,t) = (w + ecos(px)), 
p : uy(2,t) = a i{-we cos(pua) + s(n + ecos(ux))? eg + ecos(pax))° 


~ 120 
— (m+ ecos(px))}. 


etc. 


Therefore, the 3-terms Elzaki-HPM solution is: 


i 1 1 
u(a,t) = 7 + €cos(war) + aie cos(jr) + g(™ + ecos(pa))? — 156 
1 
( m+ecos(pa))’ — (4 + € cos(px))} + Teri Ga cos? (ux) + 9e® cos (tu) 
+ (36? — 32)e” cos’ (ua) + (8403 — 224m)e® cos® (ux) 
7 14 8 
4 720(u2 4 4 24 5 aad 
720(u 40" ~ 45" ip’ cos? (42) 
21 4 7 
+ 2400n( Tram pe Te cos* (11) (a) 
7 7 
— 480(€7 pu? + (—6r? + 12)? — 8x2 + an rh + e)e? cos? (ya) 
te 1 8 
— 1440r(e?p? + (—1? + 6)? 4 ce Tl 3" + €)e” cos? (ux) 
-—1 
— 1440( p76? (a? — 2) — 2u* + {=-0" — 44207} y24 4 in 
1 4 
An? TT 3m 2)e cos(ux) — 4800 (1? — 6) pe? 


384n° + 28800 — 192003 + 1° — 3207 ]t*. 


The behavior of the solution (39), by 4-terms of EHPM and its HPM solution 
are shown in Figure 3. 


Example 5.3 Considering the following Klein-Gorden equation with the 
given initial conditions: 


Utt — Use =U, u(x,0) =1+sinz, u(x,0) =0. (41) 


Taking Elzaki transform of (41), subjected to the initial condition, one obtain: 


E{(u(z,t)] = (1+ sinaz)v? + v? E[u] + Ell, (42) 


The inverse Elzaki transform implies that: 
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-60 -40 -20 0 20 40 60 


==" * EHPM ——— HPM 


Figure 3: Comparison of the results by EHPM with the HPM solution for 
example 5.2. 


u(a,t) =1+sing + BE {v?E[u] + 5 pl}. (43) 


Now, applying the homotopy perturbation method, we get: 


u(x,t) =1+sine + p{E{{v?E[u] +0? Blu}. (44) 


By taking u(x,t) = S072 p'ui(x,t) and comparing the coefficients of the 
same powers of p, result in: 


p° : uo(x,t) =1+sinz, 
t2 
p' bua (at) = 9? 
t4 
2 
: i= 
p U2(2, ) 24’ 
8 


PP :us(0,t) = zp. 


etc. 


Therefore, the 4-terms approximate series solution is: 


u(a,t) =1+sine 4 (45) 
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and this will, in the limit of infinitely many terms, yield the closed form 
solution [6], 


u(a#,t) = sina + cosht. (46) 


Example 5.4 Consider the following Klein-Gorden equation with the given 
initial conditions: 


Utt — Ure = —u, u(x,0) =0, uzt(x, 0) = a. (47) 


Taking Elzaki transform of (47), subjected to the initial conditions, one ob- 
tain: 


d2 
El(u(a2, t)] = ev? + v5 Blu] —v*Elul. (48) 
The inverse Elzaki transform implies that: 
1,27 2 
u(a,t) = ta + EB {vu — Elu] — vu" E[ul}. (49) 
Now applying the homotopy perturbation method, we get: 
-1,,2 0 2 
u(a,t) =ta+p{E "fv qe Ell —v° Elul}}. (50) 


By taking u(x,t) = 0729 p'ui(z,t) and comparing the coefficients of the 
same powers of p, one obtain: 


—at? 

Dp u(x,t) = 3! ’ 

+at® 

p u2(x, )= 5! ry 

—at" 

p? ug(a,t) = 7 ’ 
etc. 


u(x,t) = ta Al a (51) 


and this will, in the limit of infinitely many terms, yield the closed form 
solution [22], 
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u(a,t) = xsint. (52) 


6 Conclusion: 


In this paper, the Elzaki Homotopy-Perturbation method (EHPM) has been 
successfully employed to obtain the approximate analytical solutions of the 
Sine-Gorden and the Klein-Gorden equations. In example 5.3, the obtained 
result by this method is almost accurate and very near to the exact solution, 
also in example 5.4, it is observed that the (EHPM) solution yields the exact 
solution in only few iterations. Therefore, this novel iterative method has 
a bright aspect in future to obtain the approximate analytical solutions of 
ordinary and partial differential equations. 
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